Probing the Equation of State of Nuclear Matter via Neutron Star Asteroseismology 
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We general relativistically calculate the frequency of fundamental torsional oscillations of neutron 
star crusts, where we focus on the crystalline properties obtained from macroscopic nuclear models in 
a way depending on the equation of state of nuclear matter. We find that the calculated frequency 
is sensitive to the density dependence of the symmetry energy, but almost independent of the 
incompressibility of symmetric nuclear matter. By identifying the lowest-frequency quasi-periodic 
oscillation in giant flares observed from soft gamma-ray repeaters as the fundamental torsional mode 
and allowing for the dependence of the calculated frequency on stellar models, we provide a lower 
limit of the density derivative of the symmetry energy as L ~ 50 MeV. 

PACS numbers: 04.40.Dg, 21.65.Ef, 26.60.Gj, 97.10. Sj 



Observations of the global oscillations of stars play a 
significant role in probing the properties of stellar mat- 
ter, as in the case of the Sun. Studies in this direction are 
often referred to as asteroseismology. Oscillations of neu- 
tron stars are expected to give an insight to the properties 
of matter under conditions of ultra-high density. Recent 
observations of the quasi-periodic oscillations (QPOs) in 
giant flares from soft gamma-ray repeaters (SGRs) [1] 
suggest that a long-awaited evidence for the neutron star 
oscillations has been given with the help of coupling of 
gigantic magnetic fields with the solid crust. So far, three 
giant flares have been detected from SGR 0526-66, SGR 
1900-1-14, and SGR 1806-20, and, through the timing 
analysis of the X-ray afterglow, the QPOs with frequen- 
cies in the range from tens Hz up to a few kHz have been 
discovered [l|. Many theoretical attempts to explain the 
observed QPO frequencies have been done in terms of 
the torsional oscillations in the crustal region and/or the 
magnetic oscillations (e.g., Refs. If the QPOs in 

giant flares are associated with the crustal oscillations, 
the properties of inhomogeneous nuclear matter in the 
crust could be clarified lol- 1^. 

The outer part of neutron stars can be described as fol- 
lows. Below an ionic ocean in the vicinity of the surface, 
a bcc Coulomb lattice of nuclei embedded in a roughly 
uniform electron gas is considered to compose a crustal 
region. In this region, nuclei become gradually neutron- 
rich with increasing density and even drip neutrons at a 
density of about 4 x 10^^ g cm~'^. Near normal nuclear 
density, this crustal region is considered to melt into uni- 
form nuclear matter. Just before melting, roughly spher- 
ical nuclei are so closely packed that global deformations 
into rodlike nuclei could occur. As the density increases 
further, possible changes of the nuclear shape are rod, 
slab, tube, and bubble 13, These exotic nuclei are 



spherical nuclei is known to be sensitive to the empiri- 
cally uncertain density dependence of the symmetry en- 
16| . Generally, it is difficult to observationally 



15 



ergy 

probe the properties of matter in the crust, but the as- 
teroseismology could exceptionally help to constrain the 
matter properties such as nuclei present and the equa- 
tion of state (EOS) (e.g., Refs. [l^-liO]). In this Letter, 
we will show by systematic analyses that an approach 
to the QPOs in SGR giant flares in terms of the tor- 
sional shear modes in the crust could severely constrain 
the density dependence of the symmetry energy. 

We begin with the bulk energy per nucleon near the 
saturation point of symmetric nuclear matter at zero 
temperature, which can be written as a function of nu- 
cleon density n and neutron excess a 2l\ : 



W = Wo + T7-^(" 
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where wq, no, and Ko are the saturation energy, satura- 
tion density, and incompressibility of symmetric nuclear 
matter, respectively. The parameters L and So charac- 
terize the symmetry energy coefficient S{n): So = S{no) 
is the symmetry energy coefficient at rt — no, while 
L — 3no{dS / dn)n=no is the symmetry energy density 
derivative coefficient. The parameters wo, no, and So can 
be relatively easier to determine from empirical masses 
and radii of stable nuclei 22|. On the other hand, the 



often called "nuclear pasta." Not only the density region 
of the pasta phases but also the charge number of roughly 



remaining two parameters, L and Ko, are more difficult 
to determine. Here we introduce a new parameter, y, de- 
fined as y = —KoSo/i'inoL), which denotes the slope of 
the saturation line in the vicinity of a = Two of 
us (K.O. and K.I.) constructed the model for the EOS of 
nuclear matter in such a way as to reproduce Eq. ^ in 
the limit of n ^ no and a ^ 0, calculated the optimal 
density distribution of stable nuclei within a simplified 
version of the extended Thomas-Fermi theory, and ob- 
tained the values of wo, no, and So for given y (or L) 
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TABLE I: The pasta density region calculated for several sets 
of the EOS parameters. 



composed of spherical nuclei of charge Ze and number 
density n^, which can be approximately described as (27j 



y (MeV fm^) 


(MeV) L (MeV) 


m (fm-^) 


na (fm-3) 


-220 


180 


52.2 


0.060 


0.079 


-220 


230 


73.4 


0.064 


0.073 


-220 


360 


146.1 


0.066 


0.066 


-350 


180 


31.0 


0.058 


0.091 


-350 


230 


42.6 


0.063 


0.086 


-350 


360 


76.4 


0.072 


0.076 


-1800 


180 


5.7 


0.058 


0.134 


-1800 


230 


7.6 


0.058 


0.127 


-1800 


360 


12.8 


0.058 


0.118 



and ii'o by fitting the charge number, mass excess, and 
charge radius that can be calculated from the optimal 
density distribution to the empirical behavior. In a man- 
ner similar to Ref. [l^, we adopt the parameter range 
satisfying < L < 160 MeV, 180 MeV < i^o < 360 
MeV, and y < —200 MeV fm^. In fact, this parameter 
range equally well reproduces the mass and radius data 
for stable nuclei and effectively covers even extreme cases 
22{ . With such parameters, we can obtain the crust EOS 
and the equilibrium nuclear shape and size by generaliz- 
ing the Thomas-Fermi model for nuclei to matter in the 
crust. The results for the nucleon densities ni at which 
the nuclear shape changes from spherical to cylindrical 
one and n2 at which the nuclear matter becomes uni- 
form were tabulated for several sets of y (or L) and Kq 
in Table HI Note that the interval between ni and n2, 
which corresponds to the pasta region, decreases with L 
and vanishes at L ^ 100 MeV. 

We turn to the equilibrium structure of nonrotat- 
ing neutron stars, which is determined by the Tolman- 
Oppenheimer-Volkoff (TOV) equations for given EOS of 
neutron star matter. For a solution to these equations, 
the metric can be described in terms of the spherical po- 
lar coordinates r, 6, and 6 as 



ds^ ^ -e^^(-)dt^ + e2A(-)dr2 -t- r^d9^ + sin^ 9 d<t>\ (2) 

However, the EOS for the core surrounded by the crust 
is still uncertain, although there are many calculations 
based on realistic nuclear interactions and even con- 
straints from the observationally deduced neutron star 



masses and radii [231 . |24| . To avoid this uncertainty, we 
here construct the crust with the EOS models mentioned 
above by solving the TOV equations inward from the 
surface of the star for given stellar mass M and radius 
R as in Ref. 25|. We remark that the thickness of the 
crust thus constructed is consistent with a typical be- 
havior given in terms of M and R by Eq. (18) in Ref. 

126]. 

We next consider the shear modulus of a crustal part 



^ = 0.1194n,(Ze)Va, 



(3) 



where a = (3/47™^)^/'^ is the Wigner-Seitz radius. Note 
that this formula is derived in the limit of zero tempera- 
ture from Monte Carlo calculations of the shear modulus 
averaged over all directions for a perfect bcc Coulomb 
crystal of point c harg es embedded in a neutralizing uni- 
form background [28| . We will use this formula for calcu- 
lations of the torsional oscillation frequencies since effects 
of quantum zero-point motions and thermal fluctuations 
are negligible. For pasta nuclei except bubbles, the elas- 
ticity is expected to be much lower than that for spher- 
ical nuclei [29|]. Moreover, the core, which is expected 
to be composed mostly of fluids, may have a structure 
with nonvanishing elasticity J30| . Here, we simply assume 
/i = at ri > rti, as in Ref. |12l |. Within this assumption, 
the shear modulus at rt > ni is underestimated, which 
tends to lower the torsional oscillation frequencies. The 
frequencies as will be estimated below should thus be re- 
garded as lower limits basically, but still will be shown 
to play a role in constraining L. 

Since the torsional oscillations on a spherically sym- 
metric star are incompressible, the star is free from 
deformation and density variation during such oscilla- 
tions. One can thus determine the frequencies of the 
torsional oscillations with satisfactory accuracy even if 
one neglects the resulting metric perturbations by set- 
ting 5g^i, = 0, which is known as the relativistic Cowl- 
ing approximation. Within this approximation, the tor- 
sional oscillations can be described by a single perturba- 
tion variable, i.e., the angular displacement of the stel- 
lar matter, 3^, which is related to the (/)-component of 
the perturbed 4- velocity of a matter element, (Su"^, by 
5u'>' = e-"^ dty [t, r )dePt{cos 6)1 sin 9 with the ^th order 
Legendre polynomial P^. Assuming y{t,r) ~ e'"*3^(r), 
one can obtain the perturbation equation for y{r) from 
the linearized equation of motion as [sij 



y" 



- + - A' ) + ^ 
r J 
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we 



ie + 2){i-i) 



A* 



e'^y = 0, (4) 



where e and p are the energy density and pressure, and 
the prime denotes the derivative with respect to r. To de- 
termine the eigenfrequencies, we adopt the zero-traction 
condition at n = ni and the zero-torque condition at 
the star's surface. These boundary conditions reduce to 
3^' = at n = 0, ni 0, [sij. We remark that neutron su- 
perfluidity ignored here would subtract the mass density 
of superfluid neutrons from the enthalpy density e -\- p 
and hence enhance the eigenfrequencies [32|. 

First, to see the dependence of the fundamental tor- 
sional oscillations on the EOS parameters, we calculate 
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FIG. 1: (Color online) Frequency of the 1 = 2 fundamen- 
tal torsional oscillations, 0^2, plotted as a function of L for 
M = I.4M0 and R = 12 km. The horizontal dot-dashed 
line denotes the lowest QPO frequency observed from SGR 
1806-20 Q. 



TABLE II: fitting via formula © for M = 1.4Mq, where 
is given for the deviation of the calculated 0^2 from the 



fitting formula ([S]). 



R (km) 


Co (Hz) 


ci (Hz) 


C2 (Hz) 




10 


38.95 


36.14 


11.34 


1.944 


12 


34.87 


32.33 


10.14 


1.539 


14 


31.48 


29.14 


9.123 


1.239 



the corresponding eigenfrequency for a typical stellar 
model with M = 1.4Afo and i? = 12 km. Figure [D 
shows the frequency of the £ — 2 fundamental torsional 
oscillations, 0^2, calculated for the nine sets of L and Kq 
that are tabulated in Table HI together with the lowest 
QPO frequency in SGR 1806-20 fl)]. From this figure, 
one can observe that is almost independent of Kg in 
the parameter range adopted here once the stellar model 
is fixed. Such independence can be seen for various stel- 
lar models ranging R = 10, 12, and 14 km as well as 
M = I.4M0 and 1.8Mq. We can thus focus on the L 
dependence of the calculated 0^2, which arises mainly be- 
cause the nuclear charge Z decreases with L through the 
surface property [15[, leading to decrease in the shear 
modulus ([3]) with L. Since Z depends strongly on L in 
the vicinity of n = ni, this region is likely to play a role 
in constraining L via the evaluations of 0^2- To see the L 
dependence explicitly, we derive a fitting formula for 9^2 
as 



0*2 



Co - ci 



100 MeV 



+ C2 



L 



100 McV 



(5) 



where cq, ci, and C2 are the adjustable parameters that 
depend on M and R. The values of these parameters for 
M = IAMq are shown in Table HI 

The fitting formula thus obtained is exhibited in Fig. 
[21 We find from this figure that 0*2 for fixed L can be 
determined within the accuracy of ^ 20%, if R is in the 



FIG. 2: (Color online) 0^2 given by formula ((5]) for M 
IAMq and R = 10, 12, 14 km. 
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FIG. 3: (Color onhne) Same as Fig. [2] for 10 km < i? < 14 
km and 1.4Mq <M< 1.8Mq. 



range of 10-14 km. By shifting M from I.4M0 up to 
I.8M0, we find that oh for R = 10, 12, and 14 km 
decreases only by ^ 14, 10, and 9%, respectively. For 
clarity, we plot, in Fig.|3l estimated for stellar models 
ranging 10 km < i? < 14 km and 1.4M0 < M < I.SMq. 
The results are confined within the painted region, which 
is so narrow that we can constrain L as we shall see. 

Let us assume that the observed QPOs in SGR giant 
flares arise from the torsional oscillations in neutron star 
crusts and note that among many eigenfrequencies of the 
torsional oscillations, 0*2 is the lowest. Then, 0*2 would 
become equal to or even lower than the lowest frequency 
in the observed QPOs. Consequently, one can constrain 
L as L > 50 MeV from the painted region in Fig. [31 Re- 
call that the calculated 0*2 is likely to be underestimated 
because of the simplifled treatments of the shear modu- 
lus and the enthalpy density. In particular, elasticity in 
the pasta phases ignored here would act to increase 0*2 as 
long as n2 > ni 12|. At L > 50 MeV, however, n2 —ni is 
already small (see Table H]). The resultant modiflcations 
on the L constraint are thus expected to be small. 

We proceed to show the frequency of the first overtone 
of the £ = 2 torsional oscillations, it2, calculated again 
for the nine sets of L and Kq that are tabulated in Table 
[B The results with M = IAMq and R = 10, 12, 14 km 
are plotted in Fig. [H We find from this figure that 1^2 
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FIG. 4: (Color online) Frequency of the first overtone of the 
1 = 2 torsional oscillations, it2, plotted as a function of L 
for M = 1.4Af0. The solid, broken, and dotted lines are 
the results for R = 10, 12, and 14 km, while the circles, 
diamonds, and squares are those for Ko = 180, 230, and 360 
MeV. The horizontal dot-dashed line denotes one of the high 
QPO frequencies observed from SGR 1806-20 [ll. 



depends not only on L, but on Kq and R significantly, in 
contrast to the case of 0^2- The R dependence arises be- 
cause basically, 1^2 is inversely proportional to the crust 
thickness, which in turn increases as R^ [2^. The Kq de- 
pendence comes presumably from the Kq dependence of 
ni (see Table HI, but could be modified drastically once 
elasticity in the pasta phases is allowed for. This issue 
will be addressed elsewhere Q- 

In summary, we have investigated the fundamental tor- 
sional mode and the first overtone in neutron star crusts 
for various EOS and stellar models. The identification of 
the lowest QPO frequency observed from SGR 1806-20 
as 0^2 would then allow us to constrain L as L > 50 MeV. 
At present, this constraint is fairly stringent because ex- 



perimental constraints on L have yet to converge 34 1. 
While there are earlier publications that remark the sen- 
sitivity of 0^2 and 1^2 to L T2|, the present work is the 
first to provide a lower limit of L by sufficiently accurate 
and systematic calculations through the general relativis- 
tic mode evaluations and the Thomas- Fermi treatment of 
nuclei. The present constraint already suggests that the 
pasta phases, if any in neutron stars, would occur in a 
narrow density region. Neutron superfluidity and pasta 
elasticity would make the constraint on L and the pasta 
region even more stringent. 
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